Colloids in cholesterics: Size-dependent defects and non-Stokesian microrheology 
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We simulate a colloidal particle (radius R) in a cholesteric liquid crystal (pitch p) with tangential 
order parameter alignment at the particle surface. The local defect structure evolves from a dipolar 
pair of surface defects (boojums) at small R/p to a, pair of twisted disclination lines wrapping around 
the particle at larger values. On dragging the colloid with small velocity v through the medium 
along the cholesteric helix axis (an active microrheology measurement) we find a hydrodynamic drag 
force that scales linearly with v but super linearly with R - in striking violation of Stokes' law, as 
generally used to interpret such measurements. 
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Understanding the properties of colloidal particles 
moving in complex fluids is a challenging goal of soft 
matter physics and fluid dynamics. The case of a New- 
tonian fluid is very well understood, and it is accurately 
described by Stokes' law which states that the vis- 
cous force felt by a particle moving at a constant speed v 
is linearly proportional to the solvent viscosity, r]^ to v it- 
self and to the particle radius R. This fundamental and 
elegant result is exploited routinely in active microrhe- 
ology, a modern technique which consists of dragging a 
particle in a fluid, and measuring the force-velocity rela- 
tion. From these data, and the knowledge of the probe 
size, one can obtain an estimate of a fluid viscosity (ef- 
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rectly a property of the local environment. This is often 
preferrable to a bulk rheology experiment, e.g. when ad- 
dressing the jamming of dense suspensions [3]. 

What happens when a particle, instead, moves in fluid 
with broken symmetry? Here, much less is known. Sedi- 
mentation or "falling ball" experiments showed long ago 
that the drag force in nematic liquid crystals (in which 
molecular alignment is governed by a director field n) is 
anisotropic, as one might expect for a fluid of sponta- 
neously broken rotational symmetry. More recently, pas- 
sive microrheology experiments and theories have quan- 
tified the ratio of the viscosities along and perpendicular 
to the local director field, and found these to differ by 
about a factor of 2 These results were obtained in 
the linear regime where the Ericksen number Er, quan- 
tifying the ratio between viscous and elastic effects, is 
small: the particle moves "slowly" . A series of numerical 
simulations have also addressed the case of intermediate 
to high Er, corresponding to nonlinear (in v) microrheol- 
ogy. It was found that the ratio between the viscosities 
becomes much smaller in this case, and the defect struc- 
ture changes, e.g. from a Saturn ring to a dipole [El, 

Here we consider the hydrodynamics of a colloidal par- 
ticle (radius R) inside a cholesteric liquid crystal, which 
differs from the nematic by having a helical twist to 
the director field with a pitch p. (This additional or- 
der breaks translational as well as rotational symmetry.) 



We focus on the case of tangential anchoring in which the 
director field at the colloid-fluid interface lies everywhere 
parallel to the surface of the particle. We show that the 
resulting defect structure depends crucially on the ratio 
R/p. On increasing this ratio we observe a crossover from 
a dipolar configuration, with two defect patches on oppo- 
site sides of the particle, to a twisted set of disclination 
Hues of opposite chiralities wrapping the colloid. 

This crossover could profoundly influence the many- 
body self-assembly of colloids within the liquid crystal to 
form organized structures 0: we defer this issue to fu- 
ture work. Here we address a much simpler problem that 
is nonetheless directly addressable in the laboratory. We 
model an active microrheology experiment, in which we 
drag an isolated colloid through the liquid crystal, defin- 
ing two effective viscosities for motion along and perpen- 
dicular to the cholesteric helix. We find that, in contrast 
to the nematic case, the ratio between these two effective 
viscosities, even in the small Er regime, depends strongly 
on the size of the probe. This is because the effective 
drag on a particle moving along the cholesteric helix is 
superlinear in its size R (with exponent ~ 1.7), in striking 
violation of Stokes' law. This presents an instructive case 
in which microrheology delivers a probe-size dependent 
"viscosity" , essentially unrelated to any bulk macroscopic 
value, and only indirectly related to the material param- 
eters of the medium even as defined on the length scale 
R. Finally, we address the case of intermediate Erick- 
sen number (Er ^1), where the disclinations wrapping 
the dragged colloid are displaced downstream to form a 
double twisting disclination wake. 

Methodology: The thermodynamics of the cholesteric 
solvent is determined by the Landau-de Gennes free en- 
ergy Its density / is expressed in terms of a (traceless 
and symmetric) tensorial order parameter Q [8] as 
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Here Aq is a constant (setting the energy scale), K is an 
elastic constant, go = 27r/p, and 7 is a temperature-like 
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control parameter governing proximity to the isotropic- 
to-cholesteric transition. In our notation Greek indices 
denote Cartesian components and summation over re- 
peated indices is implied; e^y^s is the permutation tensor. 

We employ a 3D hybrid Lattice Boltzmann (LB) al- 
gorithm [9] to solve the Beris-Edwards equations for the 
evolution of the Q tensor [H 



AQ = r(^ + itr(|£)i). 
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Here F is a collective rotational diffusion constant and 
Dt is the material derivative for rod-like molecules [sj. 
The term in brackets is the molecular field, H, which 
ensures that in the absence of fiow Q evolves towards a 
minimum of the free energy. The velocity field obeys the 
continuity equation and a Navier- Stokes equation with a 
stress tensor, Hg^, generalised to describe liquid crystal 
hydrodynamics [95. For the hydrodynamics, the colloid is 
represented by the standard method of bounce-back on 
links where the LB distributions which link solid to 
fiuid nodes in the underlying discretised space are used 
to impose the appropriate boundary conditions at each 
step. Tangential boundary conditions for n are imposed 
on the particle surface (n is free to rotate in the tangent 
plane). Order parameter variations create an additional 
elastic force acting on the particle, which we com- 
puted by integrating the stress tensor Iiai3 over the par- 
ticle surface [12], = J dSUaisOis^ where i>i3 is the local 
normal to the colloid surface. Furthermore, the particle 
may rotate due to elastic and hydrodynamic torques. 

The thermodynamics of chiral liquid crystals is deter- 
mined by the chirality Hi and the reduced temperature 
r, which are given in terms of previously defined quanti- 
ties as: K = y^l08Kq'^/Aoj and r = 27(1 - The 
physics of a colloidal particle moving in the liquid crystal 
is controlled by the Ericksen number 



Er = jivR/K, with 71 = 2q^ /T 
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which measures the ratio between viscous and elastic 
forces. In the expression for Er, 71 is the rotational 
viscosity of the liquid crystal, q is the degree of order- 
ing in the system (for the uniaxial case with director n, 
Qai3 = Q{^a^i3 — ^a/s/^))- The Erickscu number should 
determine which dynamical regime we are in. As we 
shall see, the presence of the cholesteric helix brings an- 
other control parameter into play: R/p. This changes 
the physics completely with respect to the nematic case. 

We give most of our results in simulation units 
To convert them into physical ones, we can assume an 
elastic constant of 6.5 pN, and a rotational viscosity of 
1 poise. (These values hold for typical materials, and a 
colloidal diameter of 1 /im.) In this way, one can show 
that the simulation units for force, time and velocity map 
onto 90 pN, 1 /is, and 0.03/im/s respectively. 

Results: We first discuss the disclination structure 
around the colloidal particle. In the case of planar an- 
choring, a colloid in a nematic host leads to two surface 




FIG. 1: Defect structure close to the particle. The arrow 
denotes the direction of the cholesteric helix. For a small par- 
ticle (a), R/p =1/4, there are two defect patches on opposite 
sides. These elongate to form twisted disclination pairs for 
particles with R/p > 1/2: (b) R/p = 1/2 and (c) R/p = 3/4. 



defects of topological charge +1 (boojums) jl^. In the 
case of a chiral nematic host, we find the behavior is de- 
termined by the ratio between the colloidal size and the 
cholesteric pitch. For R/p significantly smaller than 1/2, 
we observe two defect patches on opposite sides of the 
particles (Fig. la). On the other hand, if the particle 
size exceeds the pitch, then the minimum energy con- 
figuration becomes a twisted pair of disclination lines of 
opposite chirality, each of which wraps around the parti- 
cle (Fig. lb and Ic). These can be viewed as elongated 
boojums, in which a 3d disclination line of charge -hl/2 
lies adjacent to the surface terminating in a 2d defect of 
a charge +1/2 on the surface at each end. The length 
of these "chiral rings" is not constant, rather it increases 
with the size of the particle, and for very large colloids 
{R = p) further disclination rings appear at the top and 
bottom of the particle, leading to a more uniform cover- 
age of the surface. Unlike the transition between Saturn 
rings and dipoles in standard nematic colloids [1^, this 
conformation change does not depend on the degree of 
surface anchoring. It should thus be more easily observ- 
able experimentally, purely by changing the particle size 
R or the chiral fraction (and hence p) in a mixed ne- 
matogen. Just as for the nematic case, the two defect 
patterns which we found have sufficiently different geom- 
etry that they should mediate quite distinct two-body 
and many-body effective interactions, creating new av- 
enues for directing the self-assembly of colloids by tuning 
properties of the surrounding matrix 0. 



Having seen the importance of the control parameter 
R/p in determining the statics of colloids in cholester- 
ics, it is natural to ask if this length scale ratio leads to 
different physics (from the nematic case, R/p = 0) for 
their hydrodynamics as well. To address this question, 
we studied the dynamics of colloids of various radii in 
response to a constant pulling force, which is either par- 
allel or perpendicular to the axis of the cholesteric helix. 
First we consider the Er <C 1 regime. 

Fig. 2 shows the average velocity-force relations ob- 
served, for dragging perpendicular and parallel to the 
cholesteric helix, at various R/p ratios. We fitted the 



3 




0.0004 0.0008 0.0012 0.0016 

velocity 

FIG. 2: (a) Average velocity-force characteristics for particles 
with R/p equal to 3/4, 1/2 and 1/4 respectively. The fits 
are with the formulas shown in the text. The averages were 
calculated from the last 2.5 — 5 x 10^ simulation steps. At 
steady state, the velocity along the helix shows slight periodic 
fluctuations around the average value. 



data by means of the following formula for the drag force: 

f{v) = fy^ a{R)vR (4) 

where fy and a{R)^ with dimensions of force and viscos- 
ity respectively, are parameters of the fit. We allow for a 
finite "yield force" fy , which can be seen as the microrhe- 
ological analogue of a yield stress. Our data show tha t 
the yield force, while possibly non-zero, is very small 



so that it does not appreciably change our analysis |16 |. 
Much more striking than any yield force is the extended 
regime in which the force is linear in velocity. (This would 
be the usual criterion for accepting an active microrheol- 
ogy measurement as measuring a length-scale dependent, 
but nonetheless linear, effective viscosity.) Fig. 3 shows 
the dependence of f/v ~ aR on R^ for dragging both 
perpendicular and parallel to the cholesteric helix. We 
fitted these data with a power law, f /v ^ R^. The data 
corresponding to dragging perpendicular to the helical 
axis are comparable with those reported for nematic col- 
loids: a{R) ~ R^~^ is almost constant in agreement with 
Stokes' law [l^. The fit gives a ^ 1.07 ±0.10; (note that 
the prefactor can depend on the in-plane angle between 
force and director). In contrast, pulling along the helix 
shows very different physics, with a ~ 1.7. 

One way to discuss these results is to insist on 
Stokes' law as defining an effective microviscosity r]{R) = 
a(i?)/67r, in which case for motion along the pitch r]{R) ^ 
R^-^ . This strong size effect could be attributable to the 
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FIG. 3: \n(aR)-\n{R) relation for colloids dragged perpen- 
dicular (closed squares) and along (open squares) the helix. 
Dashed and solid lines are linear fits of \n{aR) = /3 -\-a In(i^). 
The dotted line is Stokes' law f Iv — GTrr/R, for 77 = 0.1. 



R dependence of the defect structure or to non-Stokes 
features of the flow. But unless these effects are fully un- 
derstood, active microrheology is uninformative for chiral 
nematics, in the sense that it measures a combined prop- 
erty of the probe and environment, rather than any ma- 
terial property (local or otherwise) of the medium alone. 
Remarkably, this breakdown arises at very low Er, in a 
regime where the v dependence is still linear (Fig. 2). 
Our results can be illuminated by considering the bulk 
shear rheology of cholesterics undergoing permeation flow 
along the pitch direction. In this case, in any finite sam- 
ple there is a regime of linear rheology, but the shear 
viscosity increases linearly wit h sy stem size, diverging in 
the thermodynamic limit 18|, ll9| . Crudely treating the 



colloid radius R as an effective sample size then gives 
a = 2, not far from the observed value {a ~ 1.7). 

To explore further the physics involved. Fig. 4 show di- 
rector profiles for a stationary colloid, and ones moving 
along the helical axis with Er ^ 0.03, and Er ^ 0.7. In 
the stationary case the liquid crystal accommodates the 
colloid by a local deformation involving the disclination 
pattern of Fig. Ic. Interestingly, even at very small forc- 
ing the cholesteric layers slightly bend throughout the 
whole simulation box: the elastic distortions induced by 
the colloid motion appear long ranged (Fig. 4b). Here, 
the disclination rotates around the helical axis and is 
slightly displaced in the downstream direction, but still 
remains on the particle surface. At larger Er, (Fig. 4c), 
the bending of layers is more pronounced, and the discli- 
nation is fully displaced downstream to form a double 
twisting disclination "wake" (Fig. 4d). 

In summary, we have reported lattice Boltzmann sim- 
ulations of a colloidal particle embedded in a cholesteric 
liquid crystal. Focussing on tangential anchoring of the 
liquid crystal at the colloidal surface, we have shown 
that by changing the ratio between particle size and 
cholesteric pitch it is possible to control the topology 
of the local defects imposed by the presence of the col- 
loid. For small particles, the equilibrium configuration 
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FIG. 4: Snapshots of a colloid with R/p = 3/4 in a cholesteric: 
(a) stationary; (b) moving along the helical axis (Er ^ 0.03) 
and (c) likewise with Er ^ 0.7. (d) Disclination structure 
corresponding to (c). (a)-(c) were produced with QMGA [2Q| ]. 
Colors denote the magnitude of fix (the helix is along z). 



has two defects at opposite poles of the sphere. Increas- 
ing the particle size leads to a texture with two twisted 
disclination lines of opposite polarity, wrapping around 
the colloid. Both these novel defect structures should be 
detectable in experiments with crossed polarizers. When 
pulling the particle along the cholesteric helix, we have 
found the drag force is linear in v but depends superlin- 
early on particle size, in violation of Stokes' law which 
is the usual basis for interpreting active microrheology 
experiments (2j. Laser tweezers experiments (similar to 



those in nematic [2l| and in twisted nematic cells [22|) 
should be able to test our predicted superlinear scaling. 

We hope that our results will stimulate further ex- 
periments [23I and theoretical work on the microrheol- 
ogy of colloids in cholesterics. Furthermore, similar phe- 
nomenology could be expected in various complex fluid 
phases characterised by a spatially variable order param- 
eter, e.g. smectics, blue phases, and lyotropic cubic liquid 
crystals, where one expects permeation flow to occur. 
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